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1 Introduction 



One of the main purposes of accelerators being planned and built today, is to elucidate 
the mechanism of mass generation. In the Standard Model mass is generated via an 
SU{2) Higgs field doublet, associated with the existence of a Higgs particle, whereas in 
more general models there are typically several such Higgs fields, and also more physical 
particles. 

When some candidate for the Higgs particle is discovered it becomes imperative to 
establish its properties, other than the mass. While the Standard Model Higgs boson is 
even under CP, extended models may include pseudoscalar Higgs bosons. An example of 
such a theory is the minimal supersymmetric model (MSSM) |l|], where there is a neutral 
CP-odd Higgs boson, often denoted A and sometimes referred to as a pseudoscalar. 

In the context of Higgs production via the Bjorken mechanism ||, we have recently || 
investigated how angular distributions may serve to disentangle a scalar Higgs candidate 
from a pseudoscalar one. In trying to probe the uniqueness of the scalar character of the 
Higgs boson as provided by the Standard Model, we have to confront its predictions with 
those provided by possible extensions of the Standard Model. There is also the possibility 
that CP violation may be present in the Higgs sector, as first pointed out by Weinberg Q. 
We briefly discuss some possible signals of such effects. 

Below we postulate an effective Lagrangian which contains CP violation in the Higgs 
sector. In cases considered in the literature, CP violation usually appears as a one- 
loop effect. This is due to the fact that the CP-odd coupling introduced below is a 
higher- dimensional operator and in renormalizable models these are induced only at loop 
level. Consequently we expect the effects to be small and the confirmation of presence of 
CP violation to be very difficult. Although there may be several sources of CP violation, 
including the CKM matrix ||, we will here consider a simple model where the CP violation 
is restricted to the Higgs sector and in particular to the coupling between some Higgs 
boson and the vector bosons. Specifically, by assuming that the coupling between the 
Higgs boson H and the Z has both scalar and pseudoscalar components, the most general 
coupling for the HZZ-vertex relevant for the Bjorken process may be written as || [7|] 



i 2 5 / 4 ^Gp [m| <T + f (kl, k 2 2 ) (h ■ k 2 - k$k») + V (kf, k 2 2 ) e^k lp k 2a 

with kj the vector boson momentum, j = 1, 2. The first term is the familiar CP-even 
Z^Z^ H tree-level Standard Model coupling. The second term stems from the dimension- 
5 CP-even operator Z^Z^H with Z^ v = d^Z v — d u Z p . The last term is CP odd and 
originates from the dimension-5 operator e flupo Z^Z pa H . Simultaneous presence of CP- 
even and CP-odd terms leads to CP violation, whereas presence of only the last term 
describes a pseudoscalar coupling to the vector bosons. The higher- dimensional operators 
could be radiatively induced and are likely to be small. 

Related studies on how to discriminate CP eigenstates have been reported by 
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We demonstrate that a cut on the angle between the outgoing Z and the incom- 
ing electron alters the azimuthal angular distributions significantly J|. Our study also 
addresses the problem of separating the signal from the background. Furthermore, we 
include Monte Carlo data, exploiting the background study by Barger et. al. fill, and 



demonstrate that it should be possible to verify the scalar character of the Standard Model 
Higgs after three years of running at a future linear collider. We also study the energy 
spectrum of one of the final-state fermions in the Bjorken process, as recently suggested 
12| in connection with Higgs decay via vector bosons to four fermions. We compare the 



relative usefulness of the angular and energy distributions. 



2 Distinguishing CP eigenstates 

We compare the production of a Standard-Model Higgs (h = H) with the production of 
a 'pseudoscalar' Higgs particle (h = A) via the Bjorken mechanism, 

e (Pi) e + (p,) - Z (Q) h (q 3 ) -> / (g x ) f (q 2 )h (q 3 ) . (2) 

The couplings of H and A to the vector bosons are given by retaining only the first and 
last term in ([[]), respectively. 

Let the momenta of the two final-state fermions and the initial electron (in the overall 
cm. frame) define two planes, and denote by (f> the angle between those two planes; i.e. 

(pi x Q) ■ ( qi x q 2 ) 

cos0=- — -. (3) 

Pi x Q qi x q 2 



We shall discuss the angular distribution of the cross section a, 

1 da 



a 



(4) 



both in the case of CP-even and CP-odd Higgs bosons. 

The fermion- vector couplings are given by gy and g&. As a parameterization of these, 
we define the angle x by gy = gcosx and gA = gsmx- In the present work, the only 
reference to this angle is through sin 2% (see table 1 of ref. [§]). The distributions of eq. (Q) 
take the form 



2tt dan 

a H dcj) 

a a d(p 4 



1 + at(s, m) cos0 + j3(s, m) cos 2(f), (5) 
cos 2(f). (6) 



The coefficients a(s, m) and /3(s, m) given in will at very high energies vanish as s l l 2 
and s -1 , respectively. Therefore, the Standard-Model distribution (151) will asymptotically 
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at high energies become flat, whereas the CP-odd distribution in eq. (f|) is independent 
of energy and Higgs mass. A representative set of angular distributions is given in fig. 
for the case e + e~ — > fi + fi~h for both LEP2 and higher energies, and for different Higgs 
masses. There is seen to be a clear difference between the CP-even and the CP-odd 



cases. 
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Figure 1. Angular distributions of the planes defined by incoming e~ and final-state fermi- 
ons for a CP-even Higgs particle (solid) compared with the corresponding distribution for 
a CP-odd one (dashed). Different energies and masses are considered in the CP-even 
case. We assume \fs = 200 and 500 GeV at LEP2 and NLC, respectively. 



Experimentally, however, one faces the challenge of contrasting two angular distribu- 
tions with a restricted number of events and allowing also for background. We shall here 
focus on the intermediate Higgs mass range; more specifically, we consider m < 140 GeV 
where the Higgs decays dominantly to bb. The main background will then stem from 
e + e~ — > ZZ and also e + e~ — > Z7, 77. The cleanest channel for isolating the Higgs signal 
from the background is provided by the fi + fi~ and e + e~ decay modes of the Z boson. 

Let us next limit consideration to the energy range y/s = 300 — 500 GeV, as appropriate 
for a linear collider [R], henceforth denoted NLC. We impose reasonable cuts and con- 
straints as described in ITotl; e.g. \m^+^- —rriz\ < 6 GeV and | cos9z\ < 0.6. The signal for 
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e + e _ — > ZH — > [i + fi~bb will then be larger than the background e + e _ — > — > fi + fi~bb 
by an order of magnitude. In the following we shall thus neglect the background in the 
discussion of (0) versus @. With cr(e + e~ — > ZH) ~ 200 fb and an integrated luminosity 
of 20 fb -1 a year |10[|, about 4000 Higgs particles will be produced per year, in this in- 



termediate mass range. However, following [ 10 1 we have only ~ 30 signal events e + e 



ZH — > fi + fi~bb left per year for e.g. a NLC operating at \/s = 300 GeV and a Higgs par- 
ticle of mass m = 125 GeV. In the case e + e~ — ► ZH — > e + e~bb we also have a t-channel 
background contribution from the ZZ fusion process e + e~ — > e + e~(ZZ) —>■ e + e~H . This 
contribution may be neglected at LEP energies, but it is comparable to the s-channel con- 
tribution at higher energies. However, this contribution can be suppressed by imposing a 
cut on the invariant mass of the final-state electrons, e.g. \m e + e - — mz\ < 6 GeV. Hence, 
we can effectively treat the electrons on the same footing as the muons, thereby obtaining 
a doubling of the event rate. 

Imposing the cut | cos 9z\ < b, the predictions for the azimuthal correlations of eqs. @- 
(H) get modified. For the CP-even case we find || 

a\s,m) = sin 2 X sin 2 X i — 77lk 1 \ ?n , -m 2 C(&), 



4/ £,(b)\ (s, m z , m 2 ) + 12s ml 
where X(x, y, z) = x 2 + y 2 + z 2 — 2(xy + xz + yz) is the Kallen function and 



2 ( 77 — arccos b 



((b) = ^ + VT^j, C(i) = i, (8) 

whereas for the CP-odd case 

In order to demonstrate the potential of the NLC for determining the CP of the Higgs 
particle, we show in fig. |3| the result of a Monte Carlo simulation. For this purpose we 
have used PYTHIA suitably modified to allow for the CP-odd case. The statistics 
correspond to 3 years of running^] using both the n + fi~ and e + e~ decay modes of the Z 
boson. This yields about 200 events in these channels. Although the cut b = 0.6 makes a 
increase as shown in (|7|), the cos0 term is still too small to show up in the Monte Carlo 
simulation. For y/s = 300 GeV and m H = 125 GeV, the 'bare' prediction (|) for (3 is 
0.12. The cut b = 0.6 increases it slightly to 0.14. Similarly, the '—1/4' of (||) changes 

1 The event rate is based on the Standard Model, and could be different for a non-standard Higgs 
sector. 
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Figure 2. Monte Carlo data displaying the angular distribution of events e + e~ — > Zi7 - 
l + l~bb, I = fi,e for a Standard-Model Higgs versus a CP-odd one. We have taken y/s 
300 GeV, m = 125 GeV, and an angular cut | cos 9\ < b = 0.6. 



significantly to —0.39. Consequently, the cut makes it easier to discriminate between 
the CP-even distribution and the CP-odd one. From fig. ^ we see that the individual 
angular Monte Carlo distributions are consistent with the predictions, showing that a 
three-year data sample is large enough to reproduce the azimuthal distributions. In the 
Standard-Model case the fit gives 0.92±0.07 and 0.2±0.1 for the predictions 1.00 and 0.14, 
respectively, with x 2 — 1-0. In the CP-odd case the fit gives 0.94 ± 0.07 and —0.4 ± 0.1 
for the predictions 1.00 and —0.39, respectively, with x 2 — 0-7- More importantly, since 
the cos 2(f) terms are more than 4 standard deviations away, a data sample of this size is 
sufficient to verify the scalar nature of the Standard-Model Higgs. Using likelihood ratios, 
as described in |17|], for choosing between the two hypotheses of CP even and CP odd, 
we find that less than 3 years of running suffices using similar criteria. 

An alternative test has recently been suggested by Arens et. al. [0 in the context of 
Higgs decaying via vector bosons to four fermions, where one studies the energy spectrum 
of one of the final-state fermions. Applying this idea to the Bjorken process one would 
study the energy distribution of an outgoing fermion, e.g. /i~ or e~. Introducing the 
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scaled lepton energy, x = AEi-/y/s, I = fi,e, we shall consider the energy distribution of 
the cross section with respect to this final-fermion energy, 

i £ (10) 
a ax 

both in the case of CP-even and CP-odd Higgs bosons. We are using the narrow-width 
approximation and the range of x is given by x_ < x < x + , with sx± = s + m 2 z — m 2 ± y/\. 
Here the distributions are given as second-degree polynominals in x, and, as shown in H , 
the coefficients have a non-trivial dependence on the cm. energy and the Higgs mass, 
also for the CP-odd case. A representative set of energy distributions is given in fig. 
for the case e + e~ — > fi + n~h for both LEP2 and NLC energies. There is a clear difference 
between the CP-even and the CP-odd cases. 




Figure 3. Characteristic distributions for the scaled energy of the I , I = pb,e in the 
Bjorken process e + e~ — > l + l~h. Different energies and masses are considered. 



In fig. |H we show the result of a Monte-Carlo simulation for the energy distribution 
eq. (|10D analogous to the one in fig. [| Again, a cut on the polar angle is imposed. As in 
the case of angular distributions, the cut makes it easier to discriminate between the CP- 
even distribution and the CP-odd one. Also in this case the data sample reproduces the 
predicted energy distributions. An analysis of the likelihood ratios demonstrates that less 
than 3 years of running is sufficient if we require the correct answer with a discrimination 
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Figure 4. Monte Carlo data displaying the lepton energy distribution for events e + e — * 
ZH — * l + l~bb, I = fj,,e for a Standard-Model Higgs versus a CP-odd one. We have taken 
y/s = 300 GeV and m = 125 GeV. 



by four standard deviations, but more events seem to be required than in the case of 
angular distributions. 



3 CP violation 

As previously mentioned, if we allow for both the Standard- Model and the CP-odd term 
in the Higgs- vector coupling (^), then there will be CP violation. This situation will 



be discussed here. It is similar to the case of Higgs decay discussed elsewhere fll8| . We 
discard the higher- dimensional CP-even term since it is likely to be small. 

In terms of the invariant mass Si of the fermion pair, and neglecting terms of 0((lm rj) 2 ). 
the distribution (01) can be written compactly as 



0d Sl 144^(4^)4 s 2 v v 1 v ' lj 



X (s, si, m 2 ) + 4ssi (l + 2p 2 ) + 2ss 1 p 2 cos 2(0 + 5) 
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'37T\ 2 

+ sin 2% sin 2%! ( — I y/ssi (s + s\ — m 2 ) p cos(0 + 5) 
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with a modulation function 



P = Vl + (R*^) 2 A(s,si,m 2 )/(4m|), (12) 

and an angle 



Re r)(s, si) x /X(s, si,m 2 ) 



S = arctan , — n/2 < 5 < it 12, (13) 

2mi 



describing the relative shift in the angular distribution of the two planes, due to CP 
violation. This rotation vanishes at the threshold for producing a real vector boson 
(where A = 0) and, even for a fixed value of Re r], grows with energy (because of the 
V%factor). 



This relation (|13|) can be inverted to give for the CP-odd term in the coupling: 



Re 7] = —= z tan 5. (14) 
^X(s,s h m 2 ) 

This result (|TT|) is completely analogous to the one encountered for the decay of Higgs 
particles, eq. (12) of |TBJ, if we interchange <fr and ir — <p. 



Above threshold for producing a real vector meson accompanying the Higgs particle, we 
may integrate over s\ in the narrow-width approximation. Imposing the cut | cos6z\ < b, 
the distribution of eq. M) takes the compact form 



1 + a b \s,m) p cos((p + 6) + f3 b \s,m) p 2 cos 2 (0 + 5). (15) 



2tt da 
a 

with p and S given by eqs. ( |T2"D and (|l~3D, substituting s± = m%. The details of imple- 
menting a cut in polar angle are given in 0. Any CP violation would thus show up as a 
"tilt" in the azimuthal distribution, by the amount 5. 

A representative set of angular distributions is given in fig. |3| for a broad range of 
Re 7] values. We have considered a Higgs boson of m = 200 GeV accompanied by a 
/i + /i~-pair in the final state, produced at = 500 GeV. We observe that for Re r/ < 0.1 
and Re rj > 5, the deviations from the CP-even and CP-odd distributions, respectively, 
are small. Experimentally it will be very difficult to disentangle two distributions which 
differ by such a small phase shift. Thus, observation of a small amount of CP violation 
would require a very large amount of data. This should be compared with the situation 
in fig. H and fig. ||. 

We note that the special cases rj = and \r]\ ^> 1 correspond to the CP even and CP 
odd eigenstates, respectively. Hence, the distribution ( |i~5D should be interpreted as being 
intermediate between those for the two eigenstates. 
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4 Summary 

We have addressed the problem of estimating the amount of data needed in order to 
distinguish a scalar Higgs from a pseudoscalar one at a future linear collider. This is most 
likely not possible at LEP2 due to much smaller event rates and background which is 
not easily suppressed. However, we have demonstrated that one will be able to establish 
the scalar nature of the Higgs boson at the Next Linear Collider from an analysis of 
angular or energy correlations. This particular study has been carried out for the case 
t/s = 300 GeV, m = 125 GeV. Similar results are expected in other cases as long as the 
background is small. In cases where the background can not be significantly suppressed 
a more dedicated study would be required. 

In order to establish or rule out specific models, one will also need to compare different 
branching ratios, in particular to fermionic final states. The methods proposed above 
instead deal with quite general properties of the models. 

This research has been supported by the Research Council of Norway. 
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